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This study deals with the optimisation of subcritical and supercritical laminated composite drive shafts, based on a genetic algorithm. 
The first part focuses on the modelling of a composite drive shaft. Flexural vibrations in a simply supported composite drive shaft mounted 
on viscoelastic supports, including shear effects are studied. In particular, an analytic stability criterion is developed to ensure the integrity of 
the system. The torsional strength is then computed with the maximum stress criterion, assuming the coupling effects to be null. Torsional 
buckling of thin walled composite tubes is modelled using a combination between laminate theory and Fliigge theory. In the second part, the 
genetic algorithm is developed. The last part presents a comparative study between various composite materials solutions on a helicopter tail 
rotor driveline. In particular, hybrid tubes consisting of high modulus and high resistance carbon/epoxy plies are studied. These solutions 
make it possible to replace the conventional driveline consisting of five aluminium tubes by either three subcritical composite shafts or one 
supercritical composite shaft. The saving weight is equal to 42% and 72%, respectively. This study yielded some general rules for designing 
an optimum composite shaft which can be used without any optimisation algorithms. 
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1. Introduction 

Since the 1970s, composite materials have been regarded 
as potential candidates for manufacturing drive shafts of many 
kinds because of their high specific stiffness and strength [T]. 
Previous studies on this topic have dealt mainly v^rith com- 
posite shaft design in the subcritical case, that is when the 
first critical speed is never exceeded. However, when a long 
driveline is required (in the case of helicopters, tiltrotors, tail- 
less aircrafts with twin turboprops etc.), an additional means 
of increasing drive shaft length consists in to operating above 
this first critical speed, in the so-called supercritical regime. 
The main advantage of long shafts is that they reduce the num- 
ber of bearing supports, and thus greatly decrease the main- 
tenance costs and the driveline weight. The design process is 
more complex because the shaft has to go through a critical 
speed and dynamic instability due to rotating damping can 
occur in this regime. The aeronautic applications are suitable 
for operating in the supercritical regime because the driveline 
always rotates at a nominal speed during flight. Acceleration 
and deceleration operations are carried out on the ground. 
The aim of this paper is to optimise a supercritical drive shaft 
in this practical case. 

Many different numerical optimisation methods have been 
used for designing composite drive shafts in order to reduce 
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the weight, for example. Traditional methods based on the 
gradients of continuous functions have been used by several 
authors [2l|3]|4l|5). These methods are unsuitable in the case 
of composite laminates, because many of the variables which 
have to be optimised are discrete variables (such as the num- 
ber of plies) or discontinuous (such as the ply angle in clas- 
sical manufacture process). It is therefore necessary to as- 
sume these variables to be continuous to be able to compute 
gradients. The optimisation techniques for solving problems 
involving discrete variables are known as metabeuristic meth- 
ods. For example, Gubran and Gupta (6] have used simulated 
annealing techniques based on a neighbourhood approach. 
A review of the literature shows that genetic algorithms (GA) 
(71 are well adapted to designing laminate structures. GAwere 
used for the first time to optimising a composite laminate in 
8 1 and recently in the optimisation of a flexible matrix com- 
posite drive shaft in (9J. Here, it is proposed to use a GA with 
penalisation methods to account for the constraint functions. 
Additionally, in order to reduce CPU time, all design aspects 
are handled without the use of finite element methods. 

In drive shaft applications, the choice of composite ma- 
terial is of great importance. Several authors have recom- 
mended the use of hybrid composites in the production of 
drive shafts. Brown and Rezin \T0\ and Xu et al. iTTl stud- 
ied the advantages of a mixture of glass fibres and carbon 
fibres in a modified epoxy matrix. Lee et al. ||T2] and Gub- 
ran flSl studied the design and manufacture of hybrid metal- 
lic/composite drive shafts. Here it is proposed to study a mix- 
ture of high modulus (HM) and high resistance (HR) carbon 
fibre plies. The aim is to benefit from the main advantages of 
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Nomenclature 

a, Ufj 
B, B 



bit 

c 

D,Dij 

e 
E 

Eu, E22, 
E12, E44, 
E55, Em 

f 

G 



Iy,h 



J 
K 
k 
I 

m 

mat 
N 
n 
P 

q 

r 
S 
s 
T 
u 

X, y, z 



in-plane stiffness matrix of the laminate 
and elements of the matrix (1,7 = 1,2,6) 
in-plane compliance matrix of the lamin- 
ate and elements of the matrix 
coupling stiffness matrix of the laminate 
and elements of the matrix 
discrete quantity of bits 
viscous damping 

flexural stiffness matrix of the laminate and 

elements of the matrix 

thickness 

longitudinal Young's modulus of the shaft 
longitudinal and transverse Young's mod- 
ulus, Poisson's ratio, out-of-plane (trans- 
verse / normal, longitudinal / normal) and 
in-plane shear modulus of the ply 
fitness function 

transverse shear modulus of the shaft 
constraint function 

polar area moment of inertia of the shaft 
transverse area moments of inertia of the 
shaft 

polar mass moment of inertia 

reserve factors 

stiffness 

unsupported shaft section length 
mass 

material of a ply 
discrete quantity 
number of plies 
power 

number of possible orientations in the 
staking sequence 
shaft radius 

cross-section area of the shaft 
in-plane shear strength of the ply 
axial torque 

displacement (complex or real) 
coordinates 



X, X' tensile and compressive longitudinal 

strength of the ply 

Y, Y' tensile and compressive transverse 

strength of the ply 

a orientation of the ply 

77 loss factor 

jc shear coefficient 

V Poisson's ratio of the shaft 

V distribution of the torsional modes 
p mass density of the shaft 

(p in-plane cross-section rotation 

o"ll>o"22. in-plane stress of the ply (longitudinal, 

(J12 transverse and shear) 

(B natural frequency of the torsional modes 

(1) natural frequency of the flexural modes 

n spin speed 

Subscript 

B-, B+, lower and higher backward whirl speeds 

b bearing 

buck buckling 

c critical 

dv driveline 

e external 

eq equivalent 

F-, F+ lower and higher forward whirl speeds 

f flexural modes 

G gear 

i internal 

inf inferior 

m medium 

min minimum 

nom nominal 

n number of sine modes 

ply laminate ply 

s shaft 

sup superior 

T tail rotor 

t torsional modes 

th threshold speed 

str strength 



each type of fibre. The main design considerations in the case 
of composite laminate tubes are the axial stiffness and tor- 
sional resistance. In this particular case, the plies providing 
stiffness and resistance can be considered quasi independ- 
ently (14j [TS] . HM carbon fibres which have poor mechan- 
ical properties when exposed to compression loads in partic- 
ular Ugl, can serve to maximise the axial stiffness. HR carbon 
fibres can be used to maximise the resistance to torsion loads. 

The first part of this paper presents various design aspects 
such as the rotordynamics (flexural and torsional vibrations), 
torsional strength and torsional buckling. The GA is then presen- 
ted. The last part presents a comparative study between HM 
and hybrid solutions on a helicopter tail rotor driveline presen- 



ted in the literature. 

2. Design aspects 

2. 1 . Rotordynamic analysis 

When designing supercritical shafts, it is necessary to max- 
imise the external damping in order to reduce the flexural 
imbalance responses and to increase the stability in the su- 
percritical regime. There exist several techniques for increas- 
ing this non-rotating damping in a passive way. Most tur- 
bines use hydrodynamic bearings or squeeze-film dampers, 
but they have disadvantages such as the cost, the complex- 
ity of the system and the additional instabilities introduced 
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Figure 1: A simply supported axisymmetric tubular composite shaft with rolling- element bearings mounted on viscoelastic 
supports 



; whereas rolling- element bearings do not destabilise rotors 
but provide insufficient damping. However, dissipative ma- 
terials such as elastomers have recently emerged as suitable 
materials for bearing supports 1 171 118| [19| . A corresponding 
low cost configuration, consisting of an axisymmetric com- 
posite shaft simply supported on classical rolling-element bear- 
ings mounted on viscoelastic materials, is investigated here 
(Fig.[l). 

Various approaches based on beam and shell theories have 
been used to compute the flexural critical speeds of an axisym- 
metric tubular composite shaft mounted on elastic or infin- 
itely rigid supports (20l|2T]|22l|23|. The simplest of these ap- 
proaches is called the Equivalent Modulus Beam Theory (EMBT) 
[11. It is proposed here to use EMBT in the equations of mo- 
tion previously studied in (TB) : 



' 1 -H — 1 lis" + in— lis" + — u"" + — {Us- \D.Us) = 



kG 



pS " pSl 



(1) 



where u = Uy+iUz and Mb - Uby+iubz are the cross-sectional 
displacement and the deflection of the shaft respectively (see 
the list of nomenclature for other quantities) . The above equa- 
tion contains an additional term ^ in comparison with the 
equation studied in [18] : it characterizes the effects of the 
transverse shear. 

The above equation yields the four critical speeds for the 
nth harmonic (for details of the procedure, see [18|): 
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In addition in the case of a composite shaft consisting of a 
symmetric laminate, homogenised values (EMBT) have to be 
used: £= l/aneg, G = I/fleers, v - -ai2/flii, k - 2(l-i-v)/(4-i- 
3v) andp= ms/S/ wherea = A~^. 

In the field of rotordynamics, internal damping, which is 
also referred as rotating damping, is known to cause whirl in- 
stability in the supercritical regime. In the literature, the in- 
ternal damping resulting from dissipation in the shaft mater- 
ial and dry friction between the assembled components has 
been conventionally approached using the viscous damping 
model. However most materials have a vibratory damping, 
which resembles the hysteretic damping model much more 
than viscous damping, as in the case of carbon/epoxy ma- 
terials [2?, "251. Based on a previous study using the classical 
equivalence between viscous and hysteretic damping 1 18|, the 
analytical instability criterion suitable for the shaft optimisa- 
tion can be written in the following form: 
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^0 =^ 
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and where the equivalent longitudinal loss factor denoted 
rji is computed with Adams, Bacon and Ni theory [26,^1 us- 
ing complex properties of the ply (tjh =0.11 %, 7722 =0.70% 
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Table 1: Material properties corresponding to a volume fraction of 0.6 



Material 


Abbr. 


P 


Ell 


£22 


£12 


X 


X' 


Y 


Y' 


S 


fiply 








GPa 


GPa 




MPa 


MPa 


MPa 


MPa 


MPa 


mm 


Narmco 5505 [H 


BE'' 


1965 


211 


24.1 


0.36 


1365 


1586 


45 


213 


62 


0.1321 


K63712/M10 Te' 


HM 


1700 


370 


5.4 


0.3 


1500 


470 


35 


200 


75 


0.125 


T800/G947 


HR 


1530 


162 


10 


0.3 


2 940 


1570 


60 


290 


100 


0.125 



BE : boron/epoxy. 



and Tjee = 1 • 1 %) . It suffices then to compute the lowest threshold 
speed to determine the spin speed limit of the shaft. 

2.2. Torsional vibration analysis 

Assuming that the cross sections remain plane and radii 
of these cross sections remain straight, the classical torsional 
vibration equation is: 



phip^Ghcp" 



(7) 



where (p and G are the twist angle at any cross section and 
the shear modulus, respectively. The natural torsional fre- 
quencies of the shaft are mainly due to the difference between 
the mass moment of inertia of the main gearing Jq and that 
of the tail rotor Jj. This assumption gives the boundary con- 
ditions: 



/g0(O, t)^Gh(p'{0,t) 
hip{l,t)^-Gh(P'{l, t) 



(8) 

(9) 



For the tail rotor drive shaft, it can be assumed that the 
polar mass moment of inertia of the shaft /s(= Ixpl) is small 
in comparison with that of the main gear and the tail rotor. 
The natural torsional frequencies can therefore be obtained 
directly from Eqs. j7]|9} with an equivalent method of that used 
by Lim and Darlow L3J : 




Mnel 
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[n-l)n jin-D^n^ 
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2.3. Failure strength analysis 

Only the torsional resistance of the shaft is considered. 
The transmitted torque causes only in-plane shear, which can 
be computed with classical laminate theory f28l . Contrary to 
an unsymmetrical free edge laminate plate, the tubular struc- 
ture blocks the coupling effects in the case of small displace- 
ments. This can be modelled simply by assuming the clas- 
sical coupling matrix B to be null before the inversion pro- 
cedure required to compute the strain state. 



A conservative approach in the case of helicopter drive 
shaft consists in computing only the fracture of the first ply. 
Classically, the Tsai-Wu criterion (29) can be used to account 
for differences between the tensile and compressive strength, 
which can be of great importance in the case of HM fibres (see 
Table[T): 
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(12) 



where F12 is an interaction parameter which is taken to 
be equal to 0.5. It should be noted that the Tsai-Wu criterion 
includes the transverse fracture mechanism. In the HR ma- 
terial, the transverse failure strain is approximately equal to 
0.6% while the longitudinal one is equal to 1.8%. This type 
of fracture generally has no direct effects on the fracture of 
the laminate, and this approach seems to be too conservat- 
ive. Assuming that the structure will be safe up to the occur- 
rence of the first longitudinal or shear failure, we can obtain 
a more realistic torque resistance limit El|30|. This limit can 
be computed quite simply with a maximum stress criterion 



-X'<an<X ; IctizI < S 



(13) 



Comparisons between various criteria and experimental 
data are shown in Table|2] The tubes are prevent from buck- 
ling by taking a small length to diameter ratio. Tubes N°l 
and 3 shown the value of assuming presence of a null coup- 
ling mechanism in the case of unsymmetrical laminates. The 
maximum stress criterion is in better agreement with the ex- 
perimental data than the Tsai-Wu criterion. The Tsai-Wu es- 
timation is improved by taking a double value for the tensile 
transverse resistance Y. Symmetric tube N°2 is in the pure 
shear state. All the criteria are equivalent in this case. 

We previously described 1 16| the highly non-linear com- 
pressive behaviour of pitch-based fibres. The compressive 
modulus decreases dramatically which cannot be directly mod- 
elled with a linear criterion such as the above criteria. If the 
optimised solution includes +45° and/ or -45° HM fibres plies, 
it will normally be necessary to improve the strength calcula- 
tions. As we will see below, this will in fact not be necessary 
in the present case. 
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Table 2: Comparison between torque resistance calculations on various short boron/epoxy tubes 



N° 




1 


2 


3 


Stacking sequence (from inner to outer radius) 


o 


[yu,4D,-4o,yuj 


[9U,U2,9UJ 


l9U,4D,-4o,06,9tJJ 


Outer radius x lengtii 


mm 


ZD. 4 X OU.O 


Z3.4 X OU.O 


OO.O X OUO 


Experimental jl] 


N.m 


581 


132 


4689^* 


Tsai-Wu criterion 


N.m 


167 (-71%) 


130 (3%) 


1605 (-66%) 


Tsai-Wu criterion witii B = 


N.m 


313 (-46%) 


130 (3%) 


2613 (-44%) 


Tsai-Wu criterion witii Y x2 


N.m 


285 (-50%) 


130 (3%) 


1608 (-66%) 


Tsai-Wu criterion with B = and Y x2 


N.m 


496 (-14%) 


130 (3%) 


4138 (-12%) 


Maximum stress in fibre and shear directions 


N.m 


517 (-10%) 


130 (3%) 


1610 (-66%) 


Maximum stress in fibre and shear directions with B = 


N.m 


585 (2%) 


130 (3%) 


4880 (4%) 


Buckling computed with Hayashi 31 criterion 


N.m 


1049 


1547 


13016 



" Mean value of two specimen tests. 



2. 4. Torsional buckling analysis 

The finite element method is the most frequently employed 
method of computing torsional buckling. However, it is suit- 
able here a method requiring less computing time. The other 
possibility consists in solving the buckling shell theory in the 
case of orthotropic circular cylinders. Assuming that we have 
a symmetric laminate and a sufficiently long shaft, some cri- 
teria deduced from shell theory can be obtained, in particular 
that of Hayashi 1 3 1 



Tbuck - 11 



An 



42 \ 

A22 



1/4 



D 



3/4 
22 



(14) 



where is the medium radius of the shaft, and ^n, A22, 
A12 and D22 are homogenised moduli of the laminate. Shell 
theory is used here to account for the unsymmetrical lamin- 
ate. Laminate theory is included in Fliigge's shell equations 
(32l l33l . Assuming that we are dealing with a long cylinder, 
we can solve the shell equations with the simplified displace- 
ment field proposed by Fliigge |32|. In mathematical terms, 
the torque is obtained numerically by minimising a 3 x 3 de- 
terminant. These two methods of modelling and the finite 
element method for shells are compared on unsymmetrical 
stacking sequence in Table |3] All the tubes have the same 
size and the laminate all have the same thickness. Buckling 
torque is computed in the positive and then in the negat- 
ive direction. The results obtained with the finite element 
method using ABAQUS (s4 elements) (34), which were pre- 
viously validated based on experimental results obtained by 
Bauchau et al. |35| in |25|, are taken as reference values. The 
table shows that the Hayashi criterion overestimates the buck- 
ling torque, especially in the largest unsymmetrical laminates 
(N°9-12) by up to 40%. The results obtained with Fliigge the- 
ory on long tubes with laminate theory show good agreement 
with finite element calculation giving a conservative estim- 
ate on the whole. The largest errors amount to only 8% and 
the mean error is 4%. Note that the long tube assumption in- 
vestigated in [25J, is useful when the length-to-diameter ra- 
tio is greater than 100 and the diameter-to-thickness ratio is 
greater than 40. 



2.5. Driveline mass 

The driveline is composed of shafts and intermediate sup- 
ports, which include bearings, fittings, and supports. The in- 
termediate support mass m^, can be computed with an em- 
pirical equation from Lim and Darlow 13J 



mb = 17.1288 



^nnm 



0.69 



(15) 



where Pdv is the power transmitted with the driveline (in 
W) and Hnom is the nominal spin speed (in rev / min). The 
driveline mass can then be computed using the following ex- 
pression: 



- Nf^x + N\yX m\, with 



■ pSl 



(16) 



where is the number of shafts and A/b is the number of 
intermediate supports. 

3. Shaft optimisation using a genetic algorithm 

The principle of the GA algorithm is the same as that on 
which Darwin's theory of evolution is based. At the begin- 
ning, a population is randomly created and evaluated with 
a fitness function. This evaluation gives the adaptability of 
each individual. Starting with this information, the whole 
population can be evolved using selection, crossover and muta- 
tion operators. This process is iterated up to convergence. 

The main risk of this stochastic method is that of not ob- 
taining the optimum solution. In particular, GA may tend 
to converge on local optima and may not be able to cross 
these attracting points. Another weakness of the method is 
the large amount of fitness function calculations required. This 
means that the evaluation procedure must not be very time- 
consuming. 

3.1. Individual 

An individual in this driveline optimisation procedure con- 
sists of the medium diameter rm (which can be fixed or other- 
wise), the bearing stiffness kg (fixed or not), the nominal spin 
speed Onom and the stacking sequence with various materi- 
als, symmetric or not, as in the following example : [a™'''^ x 
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Table 3: Buckling torque of carbon/epoxy laminate tubes: comparison between methods of computation (/ = 4 m, rm = 40 mm, 
£11 = 134 GPa, £22 = 8.5 GPa, Eee = £55 = 4.6GPa, £44 = 4.0 GPa, £12 = 0.29, = 1.067 mm) 



Laminate ABAQUS Fliiggewith Hayashi 



laminate tlieory criterion 



N° 




IMesh" 


Nm 


Nm 


% 


Nm 


% 


1 


[15,-15]4 


60 


150 


210 


193 


-8 


222 


6 


2 


[-15,15]4 


60 


150 


214 


197 


-8 


222 


4 


3 


[30,-30]4 


60 


150 


263 


254 


-4 


283 


8 


4 


[-30,30]4 


60 


150 


268 


259 


-3 


283 


6 


5 


[45,-45]4 


60 


150 


385 


383 


-1 


419 


9 


6 


[-45,45)4 


60 


150 


385 


382 


-1 


419 


9 


7 


[02,45,-45,45,-45,02] 


60 


150 


230 


218 


-5 


252 


10 


8 


[02,-45,45,-45,45,02] 


60 


150 


219 


208 


-5 


252 


15 


9 


[02,45,0,-45,0,45,-45] 


30 


100 


358 


342 


-4 


420 


17 


10 


[02,-45,0,45,0,-45,45] 


30 


100 


329 


315 


-4 


420 


28 


11 


[02,45,02,-45,45,-45] 


30 


100 


355 


340 


-4 


440 


24 


12 


[02,-45,02,45,-45,45] 


30 


100 


313 


300 


-4 


440 


41 


13 


[-45,-15,15,45,15,-15,-45,45] 


60 


150 


389 


375 


-4 


493 


27 


14 


[45,15,-15,-45,-15,15,45,-45] 


60 


150 


439 


449 


2 


493 


12 


15 


[15,-15,-45,-15,15,45,15,-15] 


60 


150 


219 


206 


-6 


265 


21 


16 


[-15,15,45,15,-15,-45,-15,15] 


60 


150 


241 


226 


-6 


265 


10 



" Number of circumferential elements - number of lengthwise elements. 



ni, ...,a^'^^' X rij, ...,a^^^'' x Uq] where aj , nj and matj are 
the orientation, the quantity and the material of the ply 7, re- 
spectively. Under supercritical conditions, the bearing stiff- 
ness is a necessary optimisation variable because it appears 
in the rigid mode frequencies Eq. ^ and the stability criterion 
Eq. l[5|. The shaft length corresponds to the driveline length 
divided by the number of shafts. 

There are several possible ways of modelling the chromo- 
somes of individuals. It is proposed here to fix the number 
of possible orientations in the stacking sequence, denoted q. 
This sets the size of the chromosome for a particular optim- 
isation which simplifies the crossover operations. Individu- 
als are classically represented by an array of binary numbers. 
The orientation aj can be written with 2, 3 or 4 bits to stand 
for the sets {-45, 0, 45, 90}, {-67.5, -45, -22.5, 0, 22.5, 45, 67.5, 
90} or {-75, -60, -45, -30, 15, 0, 15, 30, 45, 60, 75, 90, 105, 
120, 135, 150} (in degree units), respectively, which corres- 
pond to realistic manufacturing orientations. The quantity 
nj is written with 2 or 3 bits corresponding to the sets {1,2,3,4} 
and {1,2,3,4,5,6,7,8}, respectively. The material matj is writ- 
ten with one bit to take the advantage of both HIVI and HR 
carbon fibres, or metal and HM carbon fibres, for example. 
Lastly, fce and rm are bounded and generally encoded with 
3 bits. For example, a shaft with the following stacking se- 
quence [45™"'!, O™'"^] (i.e. q = 2), with =52mmandnnom = 
with the bearing stiffness fixed and where aj and nj are en 
coded with 2 bits, rm with 4 bits, Onom with 3 bits and mat 
with 1 bit, is defined by the following chromosome: 



1 








1 








1 


1 





1 






ai- 


45° 


«i 


=2 


mati 


a2 


=0° 


"2 


=3 


mat2 



1 





1— ' 


1 








1— ' 


1 





r,„=52e[30,( 



=4000£ [3700,4400] 



The string length is therefore simply (bit^ -1- bit„ -1- bitmat) ^ 
q + hitk,+hitr^ + hita,^^. 

3.2. Constraints and fitness 

The mass is the optimised value generally used in driv- 
eline problems (31 ID |36] . In this paper, a part of the fitness 
function is inversely proportional to the mass of one shaft. 
The fitness function depends also on the resistance, buckling 
and dynamic constraints previously investigated (n £ [\l*): 



g3 = 



str str 



1 > with Kstr < 1 

nom 
^buck Tbuck 



Tn. 



Tn 



1>0 with iiTbuckSl 



1>0 



400' 



^smin 

^|eWimifeA>0 with ^tinf.^1 

^^nom 

g5n^— With ^:tsup„<l 



and in the subcritical case 

^fsupl^^cl 



86 



0.^ 



1 > with ^Cfsupi < 1 



(17) 
(18) 
(19) 
(20) 
(21) 

(22) 



6 



or in the supercritical case 



Table 5: Search parameters of the genetic algorithm 



gyn = 1 a With Kfinfn ^ 1 



8n= — 1^0 With ^CfsupnSl 



g9- 



^^nom 
-^th^th _ 
^nnm 



1 > with ^th s 1 



(23) 
(24) 
(25) 



where g; and K are constraint functions and reserve factors, 
respectively. Eq. |T7) corresponds to the torsional strength 
constraint, which requires that the torque computed with the 
resistance criterion multiplied by the reserve factor is smaller 
than the torque required. Eq. (Tb) is used here for the tor- 
sional buckling. The other equations are those giving the dy- 
namic constraints. Eqs. (20pT correspond to the positioning 
of the nominal spin speed between torsional modal frequen- 
cies. Regarding the bending modes, the constraints depend 
on which of the two cases applies: the subcritical or super- 
critical case. In the first case, Eq. (22) corresponds to the sub- 
critical assumption, i.e. the nominal spin speed multiplied 
by the reserve factor must be smaller than the first critical 
speed. In the supercritical case, Eqs. j23 24 correspond to 
the positioning of the nominal spin speed between the flex- 
ural critical speeds, and Eq. (25) corresponds to the stability 
constraint. 

GAs can not account directly of constraint functions. This 
problem can be overcome by using a penalisation method 
consisting in deteriorating the quality of an individual that vi- 
olate one or more constraints by decreasing the fitness func- 
tion. The fitness function for mass minimisation can be writ- 
ten in the general following form 



f- 



■^7;min(0,g;) 



(26) 



where jj are the penalisation factors. The reserve factors 
and penalisation factors are given in Table|4] 

3.3. The genetic algorithm method 

3.3.1. Initialisation 

The algorithm is initialised by randomly generating a pop- 
ulation of 300-600 individuals. The number depends on the 
size of the problem. 

3.3.2. Elitism 

After evaluating the population with the fitness function, 
the two fittest individuals, also called the elites, are selected 
and kept for the next generation. 

3.3.3. Scaling, selection and crossover 

With the progression of the GA, the fitness of aU the indi- 
viduals tends to converge on that of the fittest ones. This de- 
celerates the progress of the algorithm. In order to overcome 
this problem, scaling methods such as windowing, exponen- 
tial, linear transformation or linear normalisation methods 
can be used \7\. A windowing method is used here, whereby 



Population size 


300-600 


Ciiromosome lengtii 


24-34 


Crossover probability 


90% 


Mutation probability 


10% 


Number of generation 


150-40000 



the fitness of the lowest ranking individual is subtracted from 
the fitness of each individual. Two parents are then selec- 
ted, based on their scaled fitness values and a multi-point 
crossover operation is performed. The cutting point is selec- 
ted randomly. This operation gives two children, forming the 
next generation. Note that elites can be parents. 

3.3.4. Mutation 

The mutation operation consists in randomly modifying 
the bits of the chromosomes. The probability of the muta- 
tion must be very high to obtain a highly diverse population. 
But if the mutation process is too strong, the algorithm may 
not converge on the optimum fitness. Note that elites are not 
subject to mutation operations. After the mutation, the pro- 
cess is restarted at the elitism stage until the maximum fitness 
function is reached. 

The search parameters of the GA are given in Tablejs] 

4. Case study 

A helicopter tail rotor driveline presented by Zinberg and 
Symonds 1 1 1 is investigated with the GA. The original drive- 
line with a total length of 7.41 m, which is assumed to trans- 
mit a power of 447.4 kW, is composed of five subcritical alu- 
minium alloy tubes and four intermediate supports. |Zinberg 



and Symonds| proposed to replace the conventional driveline 
by three subcritical composite shafts consisting of boron/epoxy 
material. The properties of this composite shaft can be com- 
pared with those of the aluminium one in Table[6] Note that 
the Zinberg shaft was obtained only on the basis of physical 
considerations. 

In line with Lim and Darlow |3| who studied the same 
driveline case, the mass moment of inertia of the main gear- 
ing and tail rotor were assumed to be equal to 0.94 kgm^ and 
3.76kgm^, respectively. To take the difference between the 
connections in the metallic and composite shafts into account, 
a weight penalty of 1.5 kg per composite shaft is added here. 

4. 1 . Subcritical shaft optimisation 

Subcritical shaft optimisation was first studied with the 
GA. The case of two single-materials (boron/epoxy and HM/epoxy) 
and one hybrid case (HR/HM) was examinated (see Table [T] 
for the properties of the materials). The optimisation was 
computed six times with each material to check that conver- 
gence of the GA was reached. The results of the optimisation 
are shown in Fig. |2] and the properties of the fittest individu- 
als are summarise in Table [6] in each case. Figs. 2a 2d show 
the fitness function of the highest ranking individual during 
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Table 4: Reserve factors and penalisation factors 



^^:sai 


^buck ^tsupfi 
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Table 6: Optimised boron and carbon/epoxy tail rotor driveline under subcritical conditions in comparison with the conven- 
tional aluminium tail rotor driveline 



Conv. Zinberg Optimised 



Material 






aluminium 


BE|T| 


BE 


HM 


hm'' 


HR/HM 


Number of tubes 






5 






3 






String length 




bit 






24 


24 


30 


30 










[90,45, 




fqo Oo 


[90 -99 5-1 


[90™, 45™, 


tiiUlll illllci lU 








^o,ug. 


-45 45 


45 -452 




nHM _4rHR 


outer radius) 








90] 


90] 


45] 


22. 32, -6/. bj 


q/aHMi 
90 J 


Operating speed 




rev / min 


5 540 


4320 


3 800 


4 800 


4 600 


4 400 


1st critical speed 


wi 


rev / min 


8887 


5 697° 


4 606 


5 800 


5 695 


5344 


1st torsion mode 


mi 


rev / min 


2058° 


1292° 


1065 


1534 


1254 


635 


2nd torsion mode 


02 


rev / min 


65370° 


35318° 


36428 


64 965 


59 599 


34510 


Nominal torque 


Tnom 


Nm 


771 


989 


1 124 


891 


929 


971 


Strength torque 


Tstr 


Nm 


4 925° 


4 880° 


3149 


2268 


2267 


3 349 


Buckling torque 


-H^uck 


Nm 


3 090° 


2 671° 


2 645 


2108 


2105 


2 206 


Tube length 


/ 


m 


1.482 






2.470 






Mean tube radius 


rm 


mm 


56.3 


62.84 


56 


54 


50 


46 


Tube thickness 


Bs 


mm 


1.65 


1.321 


1.19 


1.00 


1.00 


1.00 


Tubes weight 


Nsiris 


kg (%) 


13.38 


8.16 (61) 


6.09 (46) 


4.26 (32) 


3.96 (30) 


3.57 (27) 


Supports weight 




kg 


15.42 


7.71 


9.68 


8.24 


8.48 


8.75 


Weight penalty 




kg 




4.5° 


4.5 


4.5 


4.5 


4.5 


Total weight 


mdv 


kg 


28.80 


20.37° 


20.27 


17 


16.95 


16.82 


Weight saving 




kg (%) 




8.4 (29)° 


8.5 (30) 


11.8 (41) 


11.9 (41) 


12.0 (42) 



Value computed vrith presented methods, with bita = 3 



the evolution (number of generations) in six different pop- 
ulations. The legend gives the properties of the fittest indi- 
vidual in each population at the last generation. The results 
obtained show that the five conventional shafts can be re- 
placed by three subcritical boron or carbon/ epoxy shafts. The 
algorithm did not find any subcritical solutions with only two 
shafts. We noted that operating speed was much higher than 
the first natural torsional frequency in all the solutions. 



In the boron/ epoxy case (Fig. 2a , GA yielded three dif 



ferent solutions with the same fitness after 2000 generations. 
All the solutions gave the same radius, the same operating 
speed and the same plies (three 90° plies, four 0° plies, one 
-45° ply and one 45° ply) but various stacking sequence or- 
ders. They also gave the same critical speeds and the same 
resistance. The independence of critical speed computations 
from the stacking sequence order is due to EMBT. As far as 
the resistance is concerned, this independence results from 
in-plane shear loading and the assumption that uncoupled 
tension-bending is involved (B = 0). On the other hand, tor- 
sional buckling is dependent of the stacking sequence order, 
as shown in Table[3] In particular, the circumferential flexural 
stiffness of the laminate is highly significant. This explains 
the position of the 90°plies, which are located in the inner and 
outer parts of the tube. The solution with the greatest buck- 



ling torque was selected as the best individual. The stack- 
ing sequence obtained ([90°2, 0°4, 45°, -45°, 90°]) is very sim- 
ilar to that of |Zinberg and Symonds| laminate ([90°, 45°, -45°, 



0°6, 90°]), only two 0° plies were replaced by a 90° ply The 
decrease in the shaft thickness and shaft radius explain the 
slight increase in weight saving from 29% to 30% obtained 
in comparison with the conventional aluminium shaft (see 
Tableje}. The computing time required for one evolution was 
approximately equal to 50 minutes using IMATLAB [37J on a 
Xeon E5540. 

The second material tested was HM carbon/ epoxy (Fig. 2b 
Convergence was reached after 200 generations, but five dif- 
ferent solutions were obtained with the same fitness. Only 
the order between 0°, 45° and -45° plies and the operating 
speed were different. The optimum shaft stacking sequence 
maximising the resistance and buckling margins was [90°, 0°3, 
45°, -45°2, 45°]. This gave the minimum thickness authorized 
(1 mm). Due to the high level of stiffness in comparison with 
boron/ epoxy, only one 90° ply was necessary to prevent buck- 
ling and three 0° plies were required to avoid reaching the first 
critical speed. The number of +45° plies was doubled due to 
the low resistance of HIVI carbon/epoxy. The weight saving 
increased considerably in comparison with the previous ex- 
ample, reaching 41% due to several combined effects: the de- 
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-[902,45,02,02,-45,90] 
'■[902,45,02,02,-45,90] 
-[902,45,02,02,-45,90] 
[902,02,45,02, -45,90]' 
'■[902,02,45,02,-45,90] 
- [902,02,02, -45,45,90] 



rris = 2.03 kg, = 56 mm. 

n„om = 3800 rpm, Lji/Qnom = 1.21, 

T,tr /T„o,r. = 3.30, Tbuck/Tnom = 2.27 



ms = 2.03 kg, Trr — 56mm, 

n„om = 3800 rpm, i^'i/finom = 1.21, 

Tstr /Tnou, = 3.30, Tbuck/T„„„, = 2.29 



ms — 2.03 kg, rm = 56 mm, 

n„om = 3800 rpm, wi/n„om = 1.21, 

Ts„ /T„„m = 3.30, Tb„ck/Tnom = 2.35 



200 400 600 800 1000 1200 1400 1600 1800 2000 
Generation 

(a) BE (bita = 2, bitmat = 0, /m € [0.05, 0.064] mm) 



0.65 




0.6 

0.55 

? 0.5 
J 

; 0.45 

1 

J 

1 

0.4 
0.35 

0.3 
0.25 

0.2 



[90, 0, O2 , -45, 452 , -45] 

[90,02,0,-45,452,-45] 

[90,45,02,0,-452,45]' 

[90,0,02,45,-452,45]. 

[90,0,45,02,-452,45]' 
[90,02,45,0,-452,45], 



m^ = 1.42 kg, rm = 54mm, 

finom = 4800 rpm, wi/n„om = 1.21, 

Tst, /T„„n, = 2.54, Tb„ck/T„om = 2.31 



ms = 1.42 kg, rm = 54mm, 

n„om = 4600 rpm, wi/Jinom = 1.26, 

Tslr /T„„m = 2.44, Tb„ck/Tn<„„ = 2.30 



ms = 1.42 kg, rm = 54mm, 

0„om = 4800 rpm, wi/n„om = 1.21, 

Tsir /T„om = 2.54, Tb„ck/T„o,„ = 2.37 



ms = 1.42 kg, rm — 54mm, 

0„„m = 4600 rpm, ii.'i/0„om = 1.26, 

Tstr /Tnom = 2.44. Tbuck/Tno,„ = 2.25 



^ — 1.42 kg, rm — 54 mm, 
Onom = 4800 rpm, ii,'x/n„om = 1.21, 

Tstr /Tnom = 2.54, Tbuck/Tnom = 2.34 



50 100 150 200 

Generation 

(b) HIM (bita = 2, bitmat = 0, rm e 10.05, 0.064] mm) 
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0.75 
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0.65 
^ 0.6 

ifi 
CU 

c 
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<D 
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0.3 



ms = 1.32 kg, 
rm — 50 mm, 
n„„m = 4600 rpm, 
u.'i/n„om = 1.20, 
Tstr/T„om = 2.28, 
Tbuck/Tnom = 2.50 



ms = 1.32 kg, 
rm = 50 mm, 
n„om = 4600 rpm, 
Wl/Onom = 1.20, 

Ts„/T„„m = 2.28, 

Tbuck/Tnom = 2.43 



[68, -68, -23, -23, 0, 23. 23, 90] 

[90. 0, -23. 23. 23, -23, 58, -68] 

[90,23,0.-23.-23.23,58,-68] _ 

[90, -23. -23, 23, 23, 23, -23, 90] 

[90, -23, -23, 23, 23, 23, -23, 90] 

[90, -23, -23, 23, -23, 23, 23, -' 




ms = 1.32 kg, 
rm = 50 mm, 
Onom = 4600 rpm, 

Uil/finom = 1.25, 
Tstr/Tnom = 2.28, 
Tbuck/Tnom = 2.38 

ms = 1.32 kg, 
rm = 50 mm, 
Hnom = 4600 rpm, 

t^l/finom = 1.24, 
Tstr/Tnom = 2.44, 
Tbuck/Tnom = 2.27 



1000 2000 3000 4000 5000 

Generation 

(c) HlVi (bita = 3, bitmat = 0, rm e [0.04, 0.054] mm) 



, -45"', 90'""] 
, -45"', 90*""] 



)'"",45'", 05'",0! 
)"■", 45"', 05'", O! 
)""'. 05"^, 0^"" , -45"'. O!;" , 90""]^ 
)"'", O""', O2'", 0"'", -45"', 90""']. 



ms — 1.19 kg, rm — 46 mm, 
0„om = 4400 rpm, uii/n„o„ 

= 1.21. Tstr/Tnom = 3.45, 

Tbuck/T„„m = 2.27 



ms = 1.20 kg, rm = 46 mm 
Onom = 4600 rpm, wi/n„„r 

= 1.24, Tstr/Tnom = 2.26, 
Tbuck/Tnom = 2.30 



)"'",05"',05'",45"',-455', 90"'" 
)"'",05"',45"',05"',-45l!',90"'" 



ms — 1.20 kg, rm = 46 mm 
n„om = 4600 rpm, wi/Qno, 

= 1.24, Tstr/Tnom = 2.26, 
Tbuck/Tnom = 2.41 



ms = 1.21 kg, rm = 42 mm, 

Hnom = 3800 rpm, Ii,'i/n„„rr 
= 1.24, Tstr/Tnom = 3.35, 
Tbuck/Tnom = 2.30 



I 2000 4000 6000 8000 

Generation 

(d) HM/HR (bita = 2, bitmat = 1, rm e [0.04, 0.054] mm) 



Figure 2: Evolution of the best individual fitness of several shaft populations in the case of various materials, subcritical condi- 



tions and three tubes forming the Zinberg and Symonds tail rotor driveline (there were 300 individuals in each evolution with 
bit„ = 1, bitr,^ = 3, bitn^^^ ^3,q^6, finom e [3800, 5200] rev / min, esmin = 1 mm and / =2.470 m) 



crease in the density, the mean tube radius, the thickness and 
the weight of the supports (due to the increase in the operat- 
ing speed, see Eq. (Ts)). 

The optimisation of the HlVi material was then carried out 
with bita = 3 i.e. a e{-67.5°, -45°, -22.5°, 0°, 22.5°, 45°, 67.5°, 
90°} (Fig.|2c|. The chromosome length increased from 24 to 
30. This considerably increased the search-space, and hence 
the number of generations required to obtain convergence 
and the computing time (approximately 3h/ evolution). Con- 
vergence was obtained after approximately 6000 generations. 
Four different optimum individuals were obtained with the 
same fitness. All of them contained (+ and -)22.5° plies and 
most of them contained (+ and/or -)67.5° plies. The optimum 



shaft selected from four solutions was [90°,-23°2 ,23°, -23°,23°2 ,- 
68°] . This shaft did not contain 0° plies and can be manufac- 
tured with a filament vnnding process. The weight saving in- 
creased slightly in comparison with the previous case due to 
the decrease in the mean tube radius. 



The last case tested was the hybrid HR/HIM case (Fig. 2d . 
Convergence was again obtained after approximately 6000 gen- 
erations, despite the fact that only three evolutions yielded 
the optimum individual. The optimum stacking sequence 
obtained was [90°""^, 45°"^^, 0°™, -45°™, 90°™]. These res- 
ults require some simple comments. The 90° and 0° plies con- 
sisted of HM fibres because these plies determine the stiff- 
ness problems (the dynamics and buckling). The +45° plies 
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consisted of HR fibres because these plies determine the strength 
problem. The weight reduction in comparison with the HM 
case was lower than expected. In fact, the decrease in the 
shaft weight was practically balanced by the increase in the 
weight of the supports. 

The optimisation procedure was also carried out in the 
case of HR material (results not presented here). In the con- 
figuration studied here, HR material gave a fitness in between 
that obtained in the case of BE and HM material, due to its 
low density. 

In addition, it is worth noting that the number of genera- 
tions required to converge on the global optimum depended 
on the size of the search-space, as well as on the basin of at- 
traction of the local and global optima. For example in the 
case of two materials with the same chromosome length, the 
number of generations required to reach convergence increased 
from 250 to 2000 (Figs. 



2a 



2b 



4.2. Supercritical shaft optimisation 

In this part, the subcritical condition was removed. The 
optimisation was performed in the case of one single -material 
(HM) and one hybrid material (HR/HM) . The results obtained 
show that the five conventional shafts can be replaced by either 
one or two supercritical shafts (see Table|7]and Figs.|3]|4}. Con- 
trary to the subcritical optimisation, it is necessary here to 
take the first four critical speeds and the threshold speed into 
account. The support stiffness was used as a supplementary 
optimisation variable to maximise the dynamic stability mar- 
gin. Lim and Darlow |3| proposed to optimise one shaft case 
with a carbon/ epoxy material denoted here CEl. (Table[l} . These 
authors used a generalised reduced gradient method that im- 
posed continuous variables such as the fraction and the ori- 
entation of the laminate plies. Among the stacking sequences 
tested, [0°a, (^/jp, -cppny 90°y]s, the optimum one was [0°64%, 
-62°4%, 90°32%]s (Table|7}. The operating speed was above 
the third critical speed and the dynamic stability was not taken 
into account. 

In the two-tube case, the first material studied was HM 
carbon/epoxy (Fig.jSaj. Convergence was reached with three 
populations after approximately 6 000 generations. The stack- 
ing sequence of the optimum shaft (involving larger margins) 
was [90°, 45°, 0°2,-45°2, 0°, 45°]. This sequence is the same as 
in the subcritical case, only the order is slightly different. The 



shaft radius increased from 54 to 56 mm. Note that in Fig. 3a 
all the margins are particularly large (> 2). In particular, the 
operating speed was in between the first and second critical 
speeds, far above the threshold speed. Much greater weight 
saving was obtained than with the conventional aluminium 
shaft (61%) or optimum subcritical shaft. This was essentially 
due to the removal of one intermediate support. 

The second case tested was that of the hybrid HM/HR 



material (Fig. 3b . Five populations converged onto the op- 
timum fitness after 10000 generations. All the stacking se- 
quences consisted of five 0° plies in HM carbon, two 90° plies 
in HM carbon, and one (-i- or -)45° ply in HR carbon. The fit- 
test individual was [90°"^, 0°™, -45 



oHR 



>HM 



90' 



,oHMi 



This 



outcome was similar to that obtained in the subcritical case. 



only one 45°^^^ ply was replaced by one 0°^^^ ply. The weight 
saving was greater than in the HM case and reaching 63%. 

In the one-tube case, only hybrid HM/HR material was 
studied (Fig. |4}. Convergence was reached with five popu- 
lations after approximately 20 000 generations. The stacking 
sequence corresponding to the optimum individual was [90°^^*, 
0°™, -45°"'^]. It is worth noting that the number of 0° plies 
increases considerably in comparison with the two-tube case 
from 5 to 9, mainly due to the fact that the dynamic con- 
straints were met. The operating speed here was in between 
the second and third critical speeds. The AG selected HR fibres 
for the 90° ply here instead of HM fibres, which was unusual, 
possibly because the tube thickness was larger than in the 
previous cases tested, which reduced the buckling risk. The 
weight saving here amounted to 72%, which is almost equal 
to that obtained by Lim and Darlow 1 3 1 . However, this solu- 
tion is more efficient because the operating speed is above 
the second critical speed and the dynamic stability was en- 
sured. 

We carried out the above optimisation procedure with bita - 
3. After 40 000 generations with 8 populations of 600 indi- 
viduals, AG did not find a better solution. In this case, the 
computing time increased dramatically, amounting to approx- 
imately 18 hours per population. 

5. Conclusion 

In the first part of this study, various models and criteria 
are presented for designing supercritical composite drive shafts. 
In particular, the following three points are worth noting. First, 
the dynamic analysis showed the destabilising effects of in- 
ternal hysteretic damping in the supercritical region. It is there- 
fore necessary to compute this effect to ensure the integrity of 
the system at the nominal speed. Secondly, the failure strength 
of composite tubes can be computed with the maximum stress 
criterion in the longitudinal and shear directions, assuming 
coupling effects to be null. The results obtained using this 
simplified approach show good agreement with the experi- 
mental data. Lastly, the torsional buckling classical criteria 
do not take the coupling mechanism involved in unsymmet- 
rical laminates into account, and this can result in large er- 
rors in the buckling torque estimates. In this particular case, 
a laminate theory obtained with Fliigge's shell equations and 
solved with Fliigge's displacement field for long tubes gives a 
better approximation of the buckling torque. 

In the second part of this study, a GA is presented for op- 
timising applications of this kind. This method was tested on 
an example available in the literature and found to be valid. 
The computing time was quite reasonable, generally amount- 
ing to less than one day. This was possible only because ana- 
lytic results and simple methods were used. On the other 
hand, the algorithm does not ensure convergence on a single 
optimum individual but yield several solutions with practic- 
ally the same mass. This study also shows the efficiency of hy- 
brid solutions based on high modulus/high resistance carbon 
fibres. In the example presented here, the hybrid solution 
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Table 7: Optimised carbon/ epoxy tail rotor driveline under supercritical conditions in comparison with the conventional alu- 
minium tail rotor driveline 



Conv. Lim 



Material 






aluminium 




rliVi 


rln/ rlivi 


ntxl rlivi 


Number of tubes 






5 


1 


2 


2 


1 


String length 




bit 






27 


33 


34 


Stacking sequence 








[064%- 


[90,45, 


rnnHM nHM 
[90 ,0^ , 


rnnHR nHM 

[90 ,0g , 


(from inner to 








—6240/,, 


02,-452, 


A cHR nHM 
— 45 ,02 , 


— 45 ] 


outer radius) 








90320/Js 


0,45] 


90 ] 




Operating speed 


^nom 


rev / min 


5 540 


6000 


5 400 


4 800 


7 000 


1st critical speed 


oil 


rev / min 


8 887 


490 


2 696 


2 647 


1 018 


2nd critical speed 


(i>2 


rev / min 




1913*^ 


10 784 


10589 


4072 


3nd critical speed 


Mo, 


rev / min 




4 303*^ 


24 264 


23 824 


9 161 


4th critical speed 




rev / min 




7 GSO*^ 


43 136 


42 355 


16 287 


1st torsion mode 


m\ 


rev / min 




ooy 


1 ozz 


/I HQ 


A QQ 


2nd torsion mode 




rev / min 




O'tD 


A Q 

41o oZO 


1 Q 1 1 9 

lo iiZ 


OUU 


Threshold speed 


i^th 


rev / min 




_b 


Zo boo 


ZU 3 JO 




Nominal torque 


Tnom 


Nm 


771 


712 


791 


891 


610 


Strength torque 


^tsai 


Nm 


4 925'' 


1492" 


2439 


2 096 


4 352 


Buckling torque 


^buck 


Nm 


3090" 


1460" 


1963 


2137 


1657 


Tube length 


/ 


m 


1.482 


7.41 


3.705 


3.705 


7.41 


Mean tube radius 


rm 


mm 


56.3 


47.7 


56.0 


50.0 


62.0 


Tube thickness 


fis 


mm 


1.65 


1.69 


1.0 


1.0 


1.375 


Support stiffness 


fce 


kNm-1 




_b 


2 864 


2 864 


1437 


Tubes weight 


NstUs 


kg (%) 


13.38 


6.08 (45) 


4.43 (33) 


3.60 (27) 


6.65 (50) 


Supports weight 




kg 


15.42 





3.80 


4.12 





Weight penalty 




kg 




1.5 


3.0 


3.0 


1.5 


Total weight 




kg 


28.80 


7.58° 


11.23 


10.72 


8.15 


Weight saving 




kg (%) 




21.22 (74) 


17.6 (61) 


18.1 (63) 


20.7 (72) 



Value computed with presented methods. Not under consideration in the reference. 
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ms = 2.22 kg, rm = 56 mm, 

n„om = 6000 rpm, ke = 2854 kN/m. 

OnomM = 2.25. UJ2/n„om = 1.78, 
t^th/Onom = 3.89. Tsir /T„om = 2.34, 

Tb„ck/T„on, = 2.67 



ms = 2.22 kg, Tm = 56 mm, 

Onoi^ = 6000 rpm, ke = 2864 kN/m. 
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Tb„ck/Tno,„ = 2.43 
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Tbuck/Tnom = 2.34 
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Figure 3: Evolution of the best individual fitness of several shaft populations in the case of various materials, supercritical 
conditions and two tubes forming the Zinberg and Symonds tail rotor driveline (there were 300 individuals in each evolution 
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= 6.48 kg, rm — 66 mm, 
n„om = 7000 rpm, ke = 1437 kN/m, 

^/ui2 = l.6A.uJ3/n„„„ = 1.37, 

l^'th/finom = 1.26, Ts„/Tn„„ = 3.57, 
Tbuck/Tnom = 2.38 



6.65 kg, Tm — 62 mm, 
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1 — /uuurpiri, — ito' MM/r 

,/u!2 = 1.72,u)3/n„o„, = 1.31, 

1^'th/nnom = 1.95, Ts,r/T„ - ' ^ 
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ms — 6.65 kg, Tm — 62 mm, 

finom = 7000 rpm, ke = 1437 kN/m. 

finomM = 1.72,U)3/n„om = 1-31, 

i-'th/Onom = 1.95, Ts„/T„„„ = 7.13, 
Tbuck/T„„„ = 2.39 
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Figure 4: Evolution of the best individual fitness of several shaft populations in the case of HM/HR material, supercritical 



conditions and one tube forming the Zinberg and Symonds tail rotor driveline (there were 600 individuals in each evolution 
with bita = 2, bitmat = 1. bit„ = 2, bitr^^ = 3, bitfc^ = 3, bitn^„^ = 3, ^7 = 5, rm £ [0.052, 0.066] mm, nnom e [5600, 7000] rev / min, 
ke e [10*, 10^] N m"\ gsmin = 1 mm, Tje = 0.1 and / =7.410m) 



made it possible to replace the conventional driveline con- 
sisting of five aluminium tubes by three subcritical compos- 
ite shafts or a single supercritical composite shaft. The weight 
saving obtained was equal to 42% and 72%, respectively. The 
supercritical solution was found to be stable, operating above 
the second critical speed. In most of the cases studied, the al- 
gorithm confirmed the general rules, previously established 
without the use of an optimisation algorithm 1 38 1 for defin- 
ing the stacking sequence of hybrid solutions: 

1. +45° HR carbon/epoxy plies should be used in order to 
maximise the torque resistance, in variable proportions 
ranging between -i-45° and -45°, depending on maximum 
torque direction; 

2. 0° HM carbon/epoxy plies should be used in order to 
maximise the axial stiffness and minimise the axial damp- 
ing; 

3. 90° HM carbon/epoxy plies should be used far from the 
middle surface in order to maximise the torsional buck- 
ling torque; 

4. the laminate does not generally have to be symmetrical. 
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